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We numerically compute the temperature dependence of spin structure factor and thermodynamic 
quantities in the antiferromagnetic quantum Ising model in the pyrochlore lattice. This model 
exhibits spin disorder ground state with exponentially-decayed spin correlation. We reproduce the 
temperature dependence of the pinch point structure in the neutron scattering experiment and 
correct entropy obtained from the measurement of the specific heat. 



I. INTRODUCTION 

Searching for magnetic monopoles is a historical and 
long persisting interest in physics. This interest was 
recently invigoratingly aroused by the discovery of the 
magnetic-monopole-like excitation in the spin ice mate- 
rial. Magnetic ions in the spin ice material, for example 
DV2T12O7 and Ho2Ti207, form the pyrochlore lattice, a 
network of corner-shared tetrahedra shown in Fig. (JlJ. 
Due to the nature of high spin quantum number^, spin 
ice material is often regarded as a classical spin system. 
The magnetic-monopole-like excitation in spin ice mate- 
rial corresponds to a classical configuration of spin ori- 
entation. Before introducing the monopole, let us review 
the ground state property first. 

The ferromagnetic exchange interaction, crystal field, 
and the dipolar interaction play the most important roles 
determining the ground state in the spin ice material. 
The anisotropic ferromagnetic exchange and the crystal 
field restrict the spins to orientate in the local [111] di- 
rection, which is along the center of the tetrahedron to 
the spin site. The spin configuration of the ground state 
in one tetrahedron is the "2-in-2-out" structure where 
two spins point outward the tetrahedron and the other 
two point inward. This structure reveals two facts in the 
spin ice material. One is the nature of the high spin 
frustration. The other is the macroscopic ground state 
degeneracy. Since there can be 6 types of ground state 
configuration in one tetrahedron, the total ground state 
degeneracy grows in proportion to the number of the 
tetrahedra. The entropy contributed by the macroscopic 
degeneracy has been measured experimentally , which is 
precisely predicted by Pauling who tries to explain the 
entropy of ice^. 

The monopole-like excitation is a single spin flip from 
one of the classical spin configuration of ground state. 
The single spin flip breaks the rule of "2-in-2-out" in two 
adjacent tetrahedra creating a monopole-anti-monopole 
pair. One tetrahedron becomes " l-in-3-out" spin con- 
figuration corresponding the monopole and the other 
becomes "3-in-l-out" corresponding the anti-monopole. 
The monopole can hop away by flipping more spins. If 



only the exchange interaction and the crystal field is 
taken into account, hopping monopole does not cost more 
energy, because there are only two tetrahedra breaking 
the rule of "2-in-2-out" in the hopping process. In this 
case, the magnetic monopole is free after its nucleation. 
If the dipolar interaction is further considered, numerical 
study shows that the interaction between monopoles be- 
comes 1/r, where r is the bond length of the separation 
distance^. 

The monopole picture seemingly explains the first or- 
der phase transition from the kagome ice state to the 
fully polarized spin state observed in Dy2Ti207 5 . As the 
magnetic field is applied in the [111] direction, the system 
undergoes a crossover to a partially polarized spin state, 
called the kagome ice state. Along the [111] direction, 
the pyrochlore lattice can be viewed as a layered net- 
work of triangular lattice and kagome lattice alternating 
to each other. The kagome ice state still obeys the rule 
of " 2-in-2-out" , but the spin orientations in the triangu- 
lar layers are quenched. Increasing the magnetic field in 
the kagome ice state generates monopole-anti-monopole 
pairs. If it is dilute, it can be regarded as a gas phase of 
magnetic monopoles. Further increment of concentration 
of monopole may lead to a liquid-gas phase transition to 
the fully-polarized state where all tetrahedra are in the 
spin orientation of "3-in-l-out"'^. 

The monopole picture to explain the magnetic proper- 
ties of the spin ice material acquired some supports from 
the neutron scattering experiments^^. An experiment of 
muon spin resonance, based on the magnetic Wien ef- 
fect, indicates that the monopople- anti-monopole gas in 
Dy2Ti207 behaves like an electrolyte^. Recently, Dun- 
siger et al. showed negative result against the monopole 
signature in Dy2Ti207^. Using the muon spin reso- 
nance, they measured the spin relaxation in the temper- 
ature region that monopole is believed to nucleate. How- 
ever, they obtained temperature independent spin relax- 
ation, that is inconsistent with the nucleation of magnetic 
monopole in Dya^O^. 

The results by Dunsiger et al. actually shaken the 
foundation of the monopole picture in the spin ice mate- 
rial. Although the monopole story is elegant and beau- 



tiful, it needs to be scrutinized, since the experiments 
supporting the monopole picture are in fact indirect evi- 
dence. One of the theoretical supports for the monopole 
picture is the 1/r 3 spin correlation in the spin ice state. 
However, there is no direct experimental evidence sup- 
porting the existence of the state with the critical spin 
correlation at zero temperature. On the contrary, there 
is experiment showing the spin freezing below 4K in 
Dy 2 Ti 2 0ff 

In this paper, we report the calculation results of the 
spin structure factor and some thermodynamic quanti- 
ties in the antiferromagnetic quantum Ising model in the 
pyrochlore lattice. The Hamiltonian is given by the fol- 
lowing 
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where a 1 are the Pauli spin matrices. This Hamilto- 
nian is relevant to the spin ice material because the ef- 
fect of ferromagnetic exchange and the crystal field can 
be effectively represented by an antiferromagnetic Ising 
model. Without the K term, the Hamiltonian describes 
a classical spin disorder phase with 1/r 3 spin correla- 
tion. For finite K, the Hamiltonian describes a coop- 
erative paramagnetic phase with exponentially-decayed 
spin correlation which adiabatically connects to the clas- 
sical paramagnetic phase in the large K limitP^l An im- 
portant feature that Eq. ([I]) is different from the classi- 
cal Hamiltonian is that the magnetic monopole is in the 
confined phased Using a model with confined magnetic 
monopole and exponential spin correlation, we are able 
to reproduce the pinch point structure and its finite tem- 
perature properties that are observed in the neutron scat- 
tering experiment!^. The numerical results in the paper 
are calculated using the Monte Carlo simulation with the 
Trotter-Suzuki decomposition. The cluster algorithm is 
applied in the imaginary time direction. 10 6 Monte Carlo 
steps with unbiased sampling and 10 5 warm-up steps are 
used in the calculation. This paper is organized as the 
following. The finite temperature results of the pinch 
point is given in Section II. The results of some thermo- 
dynamical quantities will be given in Section III. In the 
Section IV, we conclude and the summarize the relevance 
of our calculation to the spin ice problem. 



II. SPIN STRUCTURE FACTOR 

The pinch point structure predicted theoretically^ has 
been observed in the experiment of neutron scattering^. 
It is a singularity in the spin structure factor in the spin 
flip channel. At finite temperature, the singularity is 
broaden by thermal fluctuation and becomes Lorentzian 
distribution. We compute the following spin structure 
factor in the pyrochlore lattice given by 




FIG. 1. (Color online) The cubic unit cell of the pyrochlore 
lattice. There are 16 spins containing 4 tetrahedra in a unit 
cell. 



where <> is the thermal average, a, b = 1, 2, 3, 4 denote 
the basis shown in Fig. and i,j denote the fee Bravais 
lattice and 
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where e™ is unit vector of the local the easy axis. The non 
spin-flip channel of the structure factor can be computed 
by 

(i,a)(j,b);a,b>2 

The spin-flip channel Isf{Q) is obtained by subtracting 

Insf(Q) from J(Q), I SF (Q) = I(Q)-Insf(Q). A result 
of the spin structure factor in the spin-flip channel at 
K = 0.1 and T = 0.6 in the lattice of 12 x 12 x 12 unit 
cells is shown in Fig. ([2]). Since there are 16 spins in a 
unit cell, the calculation in Fig. ^ is performed in the 
system of 27648 sites. 

Our result of the spin structure factor in the spin-flip 
channel looks very similar to the experimental results by 
Fennell et a/P. In Fig. j2j), a pinch point at (0,0,2) is 
clearly observed. Due to the thermal fluctuation, the 
pinch point is broaden described by a Lorentzian distri- 
bution. We analyze the temperature dependence of the 
width of Lorentzian distribution which can be ex- 
tracted from the computational data. The results are 
given m Fig. @. The finite size effect is carefully han- 
dled. We compute for the systems of many sizes until it 
is large enough so that the value of £ is saturated for each 
temperature. The maximum number of spin computed is 
27648 in the 12 x 12 x 12 unit cells with period boundary 
condition. Our results of the temperature dependence of 
£ is consistent with the experimental results^. £(T) ~ 
exp(B/T), where B w 0.4. 
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FIG. 2. (Color online) The spin structure factor in the spin- 
flip channel is calculated at K — 0.1 and T = 0.6 (in the unit 
of J). The calculation is performed in the system consisting 
of 12 x 12 x 12 cubic unit cells. The pinch point structure at 
(0, 0, 2) is clearly seen. 
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FIG. 3. (Color online) The temperature dependence of the 
width of the Lorentzian distribution at the pinch point. The 
calculation is performed for 5 different values of temperature 
down to T = 0.4. £(T) ~ exp(B/T), where B w 0.4, with the 
R 2 value 0.9984. 



Our results indicate that the critical spin correla- 
tion is not a necessary condition for the existence of 
the pinch point structure in the neutron scattering 
data. Our model that exhibits exponentially-decayed 
spin correlatiorii" provides the same structure with the 
accurate temperature dependence. Although £ has a di- 
mension of length, it is not the typical spin correlation 
length defined in the textbook of magnetism. As a mea- 
sure of the thermal fluctuation, the thermal broadening 
of the pinch point is similar to the one of the Bragg peak 
in the neutron scattering. Therefore, the divergence of 
the £ does not imply the divergence of the spin correla- 



tion length, i.e. a critical phase, at zero temperature. 



III. THERMODYNAMICS 
A. Magnetization 

In Fig. Q, we plot the field dependence of the magne- 
tization for K — 0.1 at two different temperature T = 0.4 
and T = 0.05. The field of strength h is applied in the 
longitudinal direction. The Hamiltonian is described by 
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For each temperature, we show the results of two sys- 
tem sizes of L x L x L unit cells for L = 5 and 6. The 
results of L = 5 and L = 6 arc almost identical at both 
temperature. Namely, thermodynamic limit is reached. 
At zero field, the magnetization is zero and the system 
is spin disorder consistent with disorder by disorder^. 
As the field is applied, the system undergoes a crossover 
to state of magnetization to — 0.5 per site. Because the 
Ising variable is 1 or —1 at each site, the state of to = 0.5 
is the one that there are three sites with a = 1 and one 
site with a = —1 in every tetrahedron. Restoring the 
spin orientations by identifying the local [111] axis, it 
can be shown that the state of to = 0.5 corresponds to 
the "2-in-2-out" kagome ice state. The crossover from 
a spin disorder phase to the kagome ice state is already 
observed in the experiment!^. 
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FIG. 4. (Color online) The field dependence of the magne- 
tization at K — 0.1. The phase diagram is characterized by 
three phases: spin disorder phase at zero field, kagome ice 
phase at intermediate field, and the fully-polarized state in 
the high field. The plateau structure is also observed in ex- 
periments. Two steps of transition between those states are 
crossover similar to the classical nearest-neighbor model for 
the spin ice material. 

The crossover can be understood by the nearest- 
neighbor model of spin ice. Both the to = and the 



4 



to = 0.5 are the phases with spin gap. The field in the 
[111] direction does not close the gap, and therefore the 
transition is a crossover. Since the kagome ice state is a 
gapped state, it is robust against the magnetic field. The 
field dependence of the magnetization exhibits a clear 
plateau structure. The width of the plateau depends on 
the temperature. The lower the temperature is, the wider 
the plateau becomes. 

As field is further applied, the system undergoes an- 
other crossover from the kagome state to the fully- 
polarized state at both T — 0.4 and T = 0.05. It is a 
state of the Ising variable a = 1 at all spin sites. Restor- 
ing the spin orientation, it corresponds to "3-in-l-out" 
spin polarized state. Our result is consistent with the 
nearest-neighbor model but cannot explain the first or- 
der phase tran sition observed in the experiment. In the 
experiments^, the nature is the transition is different at 
T = 0.4 and T = 0.05. It is a crossover at T = 0.4, and 
it is a first order phase transition at T = 0.05. The first 
order phase transition is believed to the attributed by the 
dipolar interaction^ as claimed in Ref. [5] . 

B. Magnetic susceptibility 

In Fig. dHj), we compute for the temperature depen- 
dence of the inverse of the longitudinal susceptibility. 
The calculation is taken in the system of 10 x 10 x 10 
unit cells containing 16000 spins for 6 different values of 
K . In the high temperature region, all systems have the 
same susceptibility obeying the Curie- Weiss law. Since 
the magnetic susceptibility goes to zero at zero temper- 
ature in the K = case, the inverse of the susceptibility 
deviates from the Curie- Weiss law and diverges in the 
low temperature. For large K, the magnetic suscepti- 
bility obeys the Curie- Weiss law perfectly. The suscep- 
tibility at zero temperature becomes finite for finite K. 
However, due to the numerical precision, this statement 
is firmly confirmed only for K > 0.05. For K = 0.01, the 
susceptibility in the low temperature is too small beyond 
the precision of the machine. Similar to the antiferro- 
magnetic quantum Ising model in the kagome latticeP, 
the smooth transformation of the magnetic susceptibil- 
ity in the variation of K is consistent with the analysis 
in Ref. [11]. Namely, the system undergoes a crossover 
from the cooperative paramagnetic state to the classical 
paramagnetic state in the large K limit. 



C. Specific heat and entropy 

In Fig. (|6^) , we show the results of temperature depen- 
dence of Cv/T for K = 0.01, 0.05, and 3. The calculation 
is performed for the temperature down to T = 0.1 in the 
systems of 4 x 4 x 4 unit cells. The finite size effect is 
carefully considered. Although the results are not shown 
here, the system of 4 x 4 x 4 unit cells is large enough 
representing the thermodynamic limit. The entropy is 
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FIG. 5. (Color online) The temperature dependence of the 
inverse of the longitudinal magnetic susceptibility for K = 
0.01, 0.05, 0.1, 1, and 3 in the system of 10 x 10 x 10 unit 
cells. 

also computed. In Fig. ^jp), we show the results of the 
temperature dependence of the entropy. For K = 0.01 
and 0.05, the entropy in high temperature saturate at 
In 2 — \ In | , where \ In | is the residual entropy of the 
spin ice state. 

For finite K, the macroscopic ground state degeneracy 
is lifted due to the quantum fluctuation. The value of K 
introduces an energy scale of lifting. As the lowest tem- 
perature to be T = 0.1 is considered, the quantum Ising 
model gives the same results as the classical modeP for 
K = 0.01 and 0.05. We note that the error of entropy 
mainly comes from the specific heat in the low temper- 
ature. Not only the specific heat has larger error in the 
low temperature, but also the variation is sharper and 
requires finer grids of temperature in the integration of 
Cv/T. 

For K — 3, the system becomes classical paramagnetic 
phase. The entropy starts from zero and saturates at 
In 2 as temperature increases. Due to the limitation of 
the algorithm, the error of specific heat is significant in 
low temperature, because the specific heat decreases to 
zero exponentially as temperature approaches to zero. 
Numerically, it causes significant error. On the other 
hand, the specific heat at K — 3 decreases to zero slowly 
in the high temperature, leading to some error in the 
calculation of entropy. 



IV. DISCUSSION AND CONCLUSION 

In this section, we discuss the relevancy of our calcu- 
lation to the physics of spin ice material. First of all, 
we stress that we are not aiming at explaining all the 
experimental discovery in the spin ice system. The dipo- 
lar interaction which plays a significant role in the real 
systems is neglected here. The purpose of this paper is 



5 





T 1 1 1 1 

A 


L=4K=0 
L=4K=3 


50 

00 !•••»—! 


- 
















! 














- 












/ Vx. •-». 
It N!b< » 


* *- 













0123456789 10 

Temperature (in unit of J) 













0.7 


ln2 








0.6 










0.5 


ln2-(ln(3/2) )/2 


















0.4 
0.3 










0.2 


1 M 








0.1 


# 




L=4K=0.01 
L=4K=0.50 
L=4K=3.00 















0123456789 10 

Temperature (in unit of J) 



FIG. 6. (Color online) (a) The temperature dependence of 
Cv/T for K = 0.01, 0.05, and 3. For the purpose of inte- 
gration, the data of many values of temperature is needed. 
Therefore, the calculation is performed in the system of 
4x4x4 unit cells for efficiency, (b) Temperature depen- 
dence of the entropy. 



to provide a model of monopole confinement to repro- 
duce experimental results that are pronounced to be the 
evidence of deconfinement of magnetic monopole. 

The existence of the pinch point structure in the spin 
structure factor in the spin flip channel is thought to 
support the existence of the critical spin ice state. How- 
ever, our results show that the thermal broadening at the 
pinch points is not different from the thermal broadening 
at the Bragg peak measured at finite temperature. The 
vanishing of the width of the Lorentzian distribution does 
not imply the divergence of the spin correlation length. 



Furthermore, in the experiment of Chang et atf^, the 
pinch point structure is shown to be independent of the 
doping of non-magnetic ions. In the Ho2_ x Y :i ;Ti207, the 
magnitude of the pinch point is normalized showing al- 
most identical structure for x = 0, 0.3, and 1. It is well 
known that the doping effect usually changes the spin 
correlation dramatically. One of the best examples is the 
high temperature superconductor. There are also nu- 
merous transition metal oxides exhibiting diverse phase 
diagram due to doping, especially when the system is in 
the criticality. There are two possible scenarios here: 1) 
If Ho2Ti207 is in the critical spin ice state, the spin cor- 
relation should change dramatically due to the doping. 
At the 50% doping, where the system is not critical, the 
pinch point structure is still robuslP^. 2) Ho2Ti 2 07 is 
not in the criticality. Doping with non-magnetic ion, the 
system is not critical, too. In both cases, the existence of 
the pinch point structure does not imply the criticality 
at zero temperature. 

Our result of the specific heat is also consistent with 
Bramwell et alW. The saturation of the entropy reaches 
In 2 — | In | at high temperature. The missing entropy is 
contributed by the macroscopic ground state degeneracy 
but lifted by the quantum fluctuation. 

In the field dependent magnetization, the longitudi- 
nal magnetic field is applied to the local [111] direction. 
Since there are four different local axes in the funda- 
mental tetrahedra, there are four local [111] directions 
different site by site in one tetrahedron. Although our 
result of field dependent does not have experimental re- 
alization, it serves as an important result comparing with 
the classical nearest-neighbor model of spin ice. 

The cooperative paramagnetic state at finite but small 
K can be adiabatically connected to the classical para- 
magnetic state in the large K limitP} Our calculation of 
the temperature dependence of the magnetic susceptibil- 
ity for various K is consistent with that result. 

Finally, our model at small K, exhibiting the co- 
operative paramagnetic phase with confined magnetic 
monopole, has the same spin structure factor observed 
in the experiments, giving another counter example for 
the claimed existence of magnetic monopole. As sug- 
gested by Dunsiger et ai, we may need another term in 
the Hamiltonian to describe real spin ice material. Our 
results suggest that the quantum fluctuation should be 
seriously considered. 
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